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We analyze the role played by Long Distance Symmetries within the context of the Similarity 
Renormalization Group (SRG) approach, which is based on phase-shift preserving continuous uni- 
tary transformations that evolve hamiltonians with a cutoff on energy differences. We find that 
there is a SRG cutoff for which almost perfect fulfillment of Wigner symmetry is found. We discuss 
the possible consequences of such finding. 



I. INTRODUCTION 

The use of effective interactions in Nuclear Physics 
is rather old (for a review see e.g. [U [5] and references 
therein) . The basic idea is to emphasize the role of the 
physically relevant degrees of freedom, which in the case 
of nucleons in finite nuclei depends very much on the 
relevant energy scale or equivalently on the shortest de 
Broglie wavelength sampling the interactions. Moreover, 
the explicit effects of the (hard) core which traditionally 
induce strong short distance correlations actually depend 
on detailed and accurate knowledge of the interaction at 
rather short distances (see e.g. Ref. [3] and references 
therein) . 

While the issue of scale dependence is best formulated 
within a renormalization framework [3] , Wilsonian meth- 
ods have only been seriously considered as an insightful 
technique in the study of nuclear forces about a decade 
ago (for a balanced review see ^ and references therein) . 

Some years ago Glazek and Wilson [51 [7] and indepen- 
dently Wegner [8] showed how high-momentum degrees 
of freedom can decouple while keeping scattering equiv- 
alence via the so-called similarity renormalization group 
(SRG) . The SRG is a renormalization method based on a 
series of continuous unitary transformations that evolve 
hamiltonians with a cutoff on energy differences. Such 
transformations are the group elements that give the 
method its name. Viewing the hamiltonian as a matrix 
in a given basis, the similarity transformations suppress 
off-diagonal matrix elements as the cutoff is lowered, forc- 
ing the hamiltonian towards a band-diagonal form. An 
important feature of the SRG is that all other operators 
are consistently evolved by the same unitary transforma- 
tion. Quite generally the transformation may involve any 
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number of particles and allows to generate well-behaved 
multi-particle renormalized interactions. In Refs. 9,10], 
the main features of the SRG formalism were illustrated 
by Perry and one of the present authors (S.S.) through 
simple examples from quantum mechanics, namely the 
Schrodinger equations for non-relativistic two-body sys- 
tems in one and two dimensions with Dirac-delta contact 
potentials. 

Recently, the SRG has been applied to evolve several 
nucleon-nucleon (NN) potentials to phase-shift equiva- 
lent softer forms |llj which become extremely handy for 
many-body calculations in Nuclear Physics (for a review 
see e.g. 02] )• The role played by bound-states onto the 
SRG evolution has also been analyzed [13], and it was 
found that the SRG generator and evolution SRG scale 
may be suitably chosen to avoid that undesirable singu- 
larities sneak in ruining the low-energy properties of the 
SRG-evolved interaction. This has been recently applied 
to prevent the spurious bound-states which normally ap- 
pear in meson exchange potentials |14| . The complicated 
neutron-nucleus scattering becomes far simpler within a 
SRG perspective [TS] . The general discussion of operator 
evolution and in particular deuteron properties has also 
been carried out where a link to the Operator Product 
Expansion has been envisaged [16] . 

A great advantage of the SRG method over other ap- 
proaches such as the Mowk [IJ or the Unitary Correlation 
Operator Method (UCOM) method [H], is the straight- 
forward application to the scale dependence of few-body 
forces [m [Hj which consistently treat two-body induced 
as well as initially introduced few-body interactions (see 
Ref. [20J for a discussion on simple one-dimensional mod- 
els) . A further main advantage of SRG is the tremendous 
reduction of the many-body problem, since effectively the 
two-body interaction becomes almost diagonal and con- 
sequently the corresponding phase space gets enormously 
reduced. This of course is done at the cost of allowing 
three or even four-body forces which precisely due to the 
high-energy decoupling of the SRG remain shorter range 
than the two-body effective interaction but with a dif- 
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ferent and scale dependent strength. The connection to 
the UCOM configuration space method was discussed in 
Ref. [3T1 [55] where it was shown how a suitable rescaling 
of the radial coordinate is in fact equivalent to choose a 
static generator. The relation of the SRG to the Viowk 
method was discussed in Ref. [23] as a block diagonal 
cutoff representation with a Hilbert space with high- and 
low-momentum components where the mixing between 
off-diagonal elements becomes negligible along the SRG 
evolution. 

The momentum-space Mowk approach [53] (see e.g. 
Ref. 12 for a review) takes a Wilsonian point of view 
of integrating out high-energy components. This allows 
to obtain a self-adjoint Mowk potential from a bare poten- 
tial, V. Given a symmetry group with a generic generator 
X, a standard symmetry means that [T^, ^] = implies 
[^lowk,-^] — 0- The reverse, however, is not true. We 
define a long distance symmetry as a symmetry of the 
effective interaction, i.e. [Mowk,^] = but [V,X] ^ 0. 
From a renormalization viewpoint that corresponds to a 
symmetry of the potential broken only by short distance 
counterterms. This is discussed at length in Refs. [25r.27j 
and summarized in Ref. |28j . On the other hand, when 
interaction is itself evolved to low-energies the symmetry 
of the long-range part of the potential reappears explic- 
itly for the Mowk potential. 

In a recent work it has been made clear that the 
momentum-space Mowk approach ^24j displays a remark- 
able symmetry pattern which in the case of Wigner SU(4) 
spin-isospin symmetry with nucleons in the fundamental 
representation may be related unambiguously to large Nc 
dynamical features of QCD [29 as analyzed in detail in 
Refs. [25U28j . This possibility of linking a very distinct 
pattern of QCD to an observable feature of the NN inter- 
action is extremely intriguing. However, the symmetries 
are very well satisfied and a constructive point of view 
is to select what definitions of the effective interaction 
comply to the symmetries. Of course, we do not expect 
them to be perfect, but given the fact that in the Viowk 
(the diagonal elements) approach they work so well, it 
is worth testing other definitions. We want to analyze 
whether the currently used SRG complies to the symme- 
try pattern. In a previous work by one of the present 
authors (E.R.A.) [30] the running scale dependence of 
the effective interaction was carried out at low momenta 
with the sole input of low-energy scattering parameters 
such as scattering lengths (volumes) and effective ranges. 
It was found that the mixing induced by the tensor force 
was essential to achieve the Wigner symmetry condition. 
It was also found that Wigner symmetry sets in at a 
scale predicted just from the lowest threshold parame- 
ters. Also, it is clear from Ref. [25] that the difference in 
the ^Sq and '^Si phase-shifts comes from difference of the 
singlet and triplet scattering lengths while the effective 
interactions are similar in both channels. The realization 
of the Wigner symmetry does not arise from the size of 
the scattering length but rather from the long distance 
properties of the effective interactions. In Ref. [31J, the 



origin of the Wigner symmetry is unveiled in the frame- 
work of pionless effective field theory at leading order, 
but no mention is made concerning what happens when 
pion exchange interactions and the tensor coupling are 
included. 

The purpose of the present paper is to outline under 
what conditions can these long distance symmetries be 
displayed. We will show that this is indeed the case and 
moreover for the SRG cutoff of about 600MeV an ex- 
tremely accurate fulfillment of Wigner SU(4) symmetry 
is found. Taking into account that in the SRG we are 
at any rate preserving phase equivalence this result ac- 
tually suggests a representation of the interaction based 
explicitly on the symmetry. This of course will have some 
implications for Nuclear Structure and Nuclear Reactions 
which deserve further study. 

The paper is organized as follows. In Section |ll] we 
review the SRG approach and provide the definition of a 
long distance symmetry within such a framework. One 
good way to unveil a symmetry is to construct a set of 
sum rules where the symmetry is linearly broken. This is 
done in Section [IIII In Section HVl we show our numerical 
methods and the corresponding results. Finally in Sec- 
tion [V] we present our main conclusions and outlook for 
further work. In Appendix [Aj we review the meaning of 
both the Wigner and Serber symmetries within the NN 
context. In Appendix[B]we analyze a fixed-point solution 
of the SRG equations and its stability properties. 

II. SRG AND SYMMETRIES 

The formulation of the Similarity Renormalization 
Group (SRG) is well known [THl [Ml [33]- However, 
the equations are rather complicated (integro-differential 
non-linear coupled equations) and hence very little gen- 
eral properties have so far been deduced, so much insight 
is provided by numerical calculations or the study of sim- 
plified models. A textbook presentation is available [33] 
and a rigorous discussion has been carried out only re- 
cently [35 although the interesting case of unbounded op- 
erators which is the standard situation in Nuclear Physics 
has been left out. Therefore much understanding of these 
equations relies on numerical approaches and the use of 
a discretized momentum-space basis on a finite Hilbert 
space (see however [HI [3S| for some analytical models). 
The continuum limit will be analyzed in some detail in 
Appendix |B] 

The SRG makes a transformation which actually drives 
the system into a diagonal basis, suppressing exponen- 
tially off-diagonal elements, as we review here. Let us 
consider the evolution equation at the operator level in 
the SRG approach induced by the unitary transformation 

(1) 

with the anti-hermitian generator -qs — [Gg , Hs] , where 
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Gs is the SRG hermitian generator which will be specified 
shortly. If H^^o = H is the initial Hamiltonian one can 
easily see that when [Gg, Hg] = then the SRG equation 
has a stationary point. Conversely, stationary points of 
the SRG equation fulfill [[Gs,Hs],Hs] = 0. This means 
that the invariant subspaces of Hs eigenvectors are also 
invariant subspaces of Gg , whence a band-diagonal struc- 
ture follows. Due to the commutator character of the 
equations one has an infinite set of conservation laws. 
Indeed, using the cyclic property of the trace we get 



ds 



dHg 



nTr(ffrM^.,^s]]) =0. 



(2) 



In this paper we will use as initially suggested by Glazek 
and Wilson [Sl[7] the kinetic energy as the SRG generator, 
Gs = T. Using Hs = T + Vs get for the potential 
energy the evolution equation 



ds 



[[T.Vs].Vs] 



For this choice one has that for any n and s 
d 



ds 



Tr(y,") = 0. 



(3) 



(4) 



In appendix[B]we will also show that the diagonal matrix- 
elements of the standard scattering reaction R-matrix is 
a stationary point of the previous equation when the po- 
tential is diagonal. 

More generally, for an arbitrary operator O we get the 
equation 



dOs 
ds 



[[T,Vs],Os] 



(5) 



The purpose of the present paper is to outline under what 
conditions can the so-called long distance symmetries be 
analyzed. In particular, for a symmetry group generator 
X we have 



dXs 
ds 



[[T,Vs],Xs 



(6) 



A symmetry at a given scale s must commute with both 
kinetic and potential energy operators and thus fulfills, 
[Xs,T] = and [Xs,Vs] = 0. Using Jacobi's identity 

[[A, B],C] + [[C, A],B] + [[B, C], ^] - we get 



dXs 
ds 



[[Xs,T],Vs]~[[Vs,Xs],T]^0. 



(7) 



A standard symmetry corresponds to a vanishing of the 
l.h.s. whereas a long distance symmetry is a fixed- 
point of the evolution along the similarity renormaliza- 
tion group at a given point. We will see that such a 
symmetry pattern appears within the SRG in regard to 
Wigner and Serber symmetries (see Appendix [A] for a 
short description). The variable s has dimensions of 



energy"^ and it is customary to introduce the SRG cutoff 
A = s~^/'* which has dimensions of momentum. 

As we have already mentioned, the previous manip- 
ulations are fully justified in a finite Hilbert space and 
typically corresponds to a discretized momentum-space 
basis. The discretization effects are enhanced as the run- 
ning Hamiltonian is driven towards the diagonal form. 
So one expects to evolve to a scale where the high-energy 
components are decoupled from the dynamics but also 
where the discretization effects are also unimportant. A 
feature of the SRG is that much of the SRG scale evo- 
lution happens already at the beginning and slows down 
as one approaches the scales of interest in nuclear appli- 
cations. 

We expect the long distance symmetry (if at all) to 
happen at a given SRG scale, but this may depend on the 
SRG generator. For definiteness, we use here the SRG 
equations with the kinetic energy as the generator, which 
has T\p) — Ep\p) with Ep —p^/M. In momentum-space 
the equations read 



dVs{P,p) 



ds 



d^q 



-{Ep-E'^fVs{^,p) 
{Ep + Ep, - 2Eg) Vsiif, q}Vs{q,p) . (8) 



Note that if we take zero momentum states 



dVs{0,6) 



= -2 



d\ 
(27r)3 



E,{g\VsV,^\q) <0, (9) 



which means that zero- momentum matrix elements of the 
SRG evolved potentials always decrease. Note that this 
does not prevent that it becomes infinite at some finite 
value of s. The conservation of Tr (VJ^) means that the 
zero momentum strength increases at the cost of deplet- 
ing the high- and off-diagonal matrix elements. 



III. SYMMETRIES AND PARTIAL WAVES 
DECOMPOSITION 

A. Kinematical Symmetries 

The most general self-adjoint interaction in 
momentum-space which is invariant under parity, 
time-reversal, isospin and Galilean invariance has the 
following form [37', 

V{p',p) = Vc + fi • 7*2 Wc + {Vs +Ti-T2Ws)^l- CT2 

+ (VlS + n • T2 Wls) + i?2) • (p xp) 

+ {Vt + Ti ■ f2WT) Si2il^ - p) 
+ {Vq + Ti ■ f2WQ) Si2{l^ X p) 
+ {Vp+n-T2Wp)Si2{lf +P), 

where the tensor operator is defined as 



1 



CTi • g(T2 • g 



-q (Ti • CT2 



(10) 
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and Bi and Ti are the Pauli matrix spin and isospin oper- 
ators for i-th particle respectively. The subscripts refer 
to the central (C), spin-spin (S), tensor (T), spin-orbit 
(SL), quadratic spin-orbit (Q) and quadratic velocity de- 
pendent (P) components of the NN interaction, each of 
which occurs in an isoscalar (V) and an isovector (W) 
version. Here cti 2 and fi 2 are the usual spin and isospin 
operators of the two nucleons. Note that the so-defined 
condition of zero angular averaging for the tensor oper- 
ator / (P pS 12 {p) = is fulfilled and as a consequence 
our results will be stated in a rather simple form. Note 
also that the central parts can also be deduced by tracing 
the potential in the spin-isospin space after appropriate 
multiplication of the operators 1, S^, and Gia (see 
Appendix [A| 

An advantage of using the momentum-space basis is 
that the generalized Pauli principle can be incorporated 
directly into the potential which is antisymmetric under 
the exchange of the final states which in the CM reads, 



Vl',2';l,2(p',P) = -V2',1';1,2(-P',P) 



(11) 



where the indices represent a full Pauli spinor-isospinor 
state. These kinematical symmetries are preserved by 
the SRG equations. This means that if we have a start- 
ing potential V{p',p) which evolves into Vs{f ,p) then 
V{Rp' , Rp) evolves into Vs{Rp' , Rp). In particular, a po- 
tential remains rotational invariant along the SRG evolu- 
tion and the general form is maintained throughout the 
evolution. This of course includes the generalized Pauli 
principle (111. 



where Vs{p,p') is used as a brief notation for the pro- 
jected NN potential matrix elements. 



JLL'S-J] 



{p,p') = {p{LS)J-I\Vs\p'{L'S)J-I) , (14) 



in a partial-wave relative momentum-space basis, 
I q{LS)J; I ), with normalization such that 



1 = -/ dqq'\q{LS)J;I) {qiLS)J;I\. (15) 

The superscripts J, L{L'), S and / denote respectively 
the total angular momentum, the orbital angular momen- 
tum, the spin and the isospin quantum numbers of the 
NN state. For non-coupled channels {L = L' = J), the 
matrix elements Vs{p,p') are simply given by Vs{p,p') = 



(JJJS;I), , 



{p,p'). For coupled-channels {L,L' = J± 1), 



the Vs{PtP') represent 2x2 matrices of matrix-elements 
for the different combinations of L and L' . The advantage 
of using the partial-wave decomposition is that every sin- 
gle channel may be evolved independently from a given 
initial solution, Vs=o{p' ,p), to yield a unitarily equivalent 
potential 

V,{p\p) = I dqqUs{q',p'rVs=o{q',q)Us{q.P) (16) 



A further property is that for low-momenta we have the 
threshold behaviour 



B. Long distance Symmetries 

The off-shell T— matrix has the same decomposition 
as for the potential and hence contains 12 independent 
amplitudes but on-shell, because of the condition jp*! — 
IpI, one just gets 10 different amplitudes since one has the 
identity (Ti • (72 — (7i-na2-n+ai-m<72''>Ti+<^i'la2-l where h, 
m, I are three orthonormal vectors in three-dimensional 
space such as those proportional top —p,f +p and pf xp. 

Following the description in ref. [35] (see also [321110]) 
one can undertake the projection onto partial waves. The 
non-linear SRG equation can be simplified using the con- 
servation of angular momentum. If we just define for a 
given spin the standard partial wave decomposition 



{tf\V,^\p)^N 



E 

JMLL' 



with N — Att"^ /M we get an infinite set of coupled equa- 
tions with good total angular momentum. The three- 
dimensional structure contains central and non-central 
forces. Inserting this into the SRG equation we get the 
coupled-channel equations, 

dVs{p,p') 



ds 



{p'-p''rvs{p,p') 



- / dq (p2 +p'2 _ 2q^) Vs{p,q) Vs{q,p'){13) 



Vs{p\p) = Cl,l'{s)p\p')^ [1 + 0{p\p'^)\ (17) 

which means also that the matrix Cl^l'{s) is a decreas- 
ing function of the SRG cut-off, i.e. l'(s) is negative 
definite. 



C. Perturbation theory 

In the case of rotational invariance, the symmetry is 
preserved along the SRG trajectory. In the case of a 
long distance symmetry the symmetry breaking strength 
depends on the SRG scale. Let us thus consider the ex- 
pansion around the central solution, i.e. let us split 



Vnn ^Vo + Vi 



(18) 



where [L, Vq] = whereas [J, l^i] = and [L, Vi] ^ 0. 
The zeroth-order potential commutes with L, S,T and 
so the corresponding potential may be denoted as Vf^. 
The total potential commutes with the total angular mo- 
mentum J = L + S and S^. Of course the goodness of 
this separation depends on the SRG scale A, but we ex- 
pect it to work better the lower the scale. Inserting the 
decomposition (18) into the SRG equation we have the 
zeroth-order 
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ds 

whereas the first-order is 
ds 



+ 



d\ 



(2^)3 



{Ep + E.p, - 2E,) VP il/, q)Vj"'> {q, p) + VP {j/, q)VP {q, p) 



(19) 



(20) 



Clearly, because Vq is central we get that the linear com- 
binations are preserved as long as the non-central contri- 
bution remains a small correction. Generally, we do not 
expect it to be the case, and in particular for the initial 
potential. 



D. Sum rules 

Rather than analyzing the evolution of the generators 
we prefer to check the running of a set of sum rules based 
on first-order perturbation theory in non-central compo- 
nents of the potential, which helped to disentangle some 
correlations in TV TV fits [JTJ 02] and have been shown to 
work well at the level of phase-shifts and Fiowk potentials 
in Ref. [ISHUj. The main idea is to decompose the poten- 



tial into central and a non-central pieces, as in Eq. ( 18 ) 
and assume the non-central piece to be of the form 



Vi = L ■ SVls + SuVt , 



(21) 



One property fulfilled by the non-central interaction is 
that the trace in the spin-isospin space vanishes. This 
suggests a set of sum rules at the partial waves level. 
Using first-order perturbation theory [25l [26] we deduce 
the following linear combinations for triplet P-waves, 



1 



{V.p„+3Vsp^+5V.p,) 



72 
1 

12 



3V3 



"Pi 



V3 



'P2 



{2V.ip,+3V3p^-5V3p, 



for triplet D-waves 



1 o 

7 



'"Dj 



120 
1 

60 



(3^3^^ - 5V3D, + 2V30, 



(22) 
(23) 
(24) 

(25) 
(26) 



W3D,+^V3o,-lW3o,), (27) 



for triplet F-waves 

y-Fc = ^{hV3p^+W3p^+W3F,), (28) 

= -^{^y-F.-iyF^ + iyp,). m 

yp.s = -I^(20^^^^2 + 7y3F, - 27V3F,) , (30) 

and for triplet G-waves 
1 



y^Gc 



27 



[7V3G,+9V3G, + liyG,), 



V3 



77 



3n^ = — : 



2160 



5^3, 



9V3 



W3 



(31) 
(32) 



y-G.s 



360 



iW3G,-W3G,+UV3G,). (33) 



In terms of the previous definitions Serber symmetry 
reads 

= V3p^ = V3p^: ^V3Ha=... (34) 

whereas Wigner symmetry implies 

Vis, = v3s^. viD.^yoc^ y-G.^y^Gc- (35) 



IV. NUMERICAL RESULTS 

Here we provide some details on the numerical evolu- 
tion. In order to solve the Wegner's flow equation, we 
discretize the momentum-space using a gaussian grid of 
N mesh-points. This leads us to a system of iV^ non- 
linear coupled first-order differential equations. The sys- 
tem can then be solved with an implementation of a fifth- 
order Runge-Kutta algorithm with adaptive step. In this 
work we use N = 200 points. 

Since the potentials we are evolving through the SRG 
are all regular, we set an ultraviolet cutoff at A = 30 fm^^ 
which is beyond the point where regulated potentials van- 
ish. In principle, one could take a larger cutoff value 
provided the number of points in the grid are enough to 
ensure the convergence of the Runge-Kutta algorithm. If 
the cutoff is too large, one needs too many points in the 
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FIG. 1: Diagonal matrix-elements of the SRG-evolved Argonne AV18 potential 131 V(j),p) (in fm) as a function of the CIVI 
momentum (in fm~^) for the S, P, D, F and G partial-wave components for different values of the similarity cutoff A. Left 
Panel: Initial potential (A = oo). Central Panel: A = 2 fm~^. Right Panel: A = 1 fm~^. 
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grid and the number of coupled equation becomes too 
large to be solved in a reasonable time. 

We illustrate our points for the Argonne AV18 poten- 
tial [37 which fits not only the NN phase-shifts of the Ni- 
jmegen data base PWA ^46j, but also the deuteron elastic 
form factors and has been used quite often successfully for 
nuclear structure calculations The diagonal matrix- 
elements V{p,p) and fully off-diagonal matrix-elements 
V{p, 0) are depicted at Fig. [l] and Fig. [2] respectively. 
We compare the initial potentials and the SRG poten- 
tials evolved to A = 2 fm~^ and A = 1 fm^^. As we 
see, already above A = 2 fm^^ the potentials for the ^Sq 
and ^Si waves cross. This is an indication that Wigner 
symmetry works very well around that scale. This trend 
is also observed in even partial waves such as D and G, 
where the effect of SRG evolution becomes less impor- 
tant as the angular momentum is increased. On the other 
hand, odd partial waves provide hints of the Serber sym- 
metry as one sees that the potential is much larger 
than the ^Lq combination. 

A comparison between the Argonne AV18 poten- 
tial [13], the Nijmegen II potential [H] and the chiral 
N3L0 potential of Entem and Machleidt [45] is presented 
in Fig. [s] for the similarity cutoff A = 2 fm~^ and for the 
S- and P-waves. As we see, there is some degree of uni- 
versality, as one might expect since these potentials are 
phase equivalent, although the strength is distributed dif- 
ferently, particularly in the N3L0 chiral potential case. 

As we see, the potential for the ^S'l wave changes 
rather dramatically when going from A = 2 fm~^ to 
A = I fm~^ in contrast to the other waves where appar- 
ently a much more stable result is obtained. This feature 
is due to the use of the simple generator rjs = [T,Hs], 
where T is the kinetic energy. In Wegner's original for- 
mulation [S] the generator rj = [Ds,Hs] was used, where 
Ds is the diagonal part of Hg. Actually, as shown by 
Glazek and Perry [T3] the generator r/s = [T, Hg] can 
produce divergencies in theories with bound-states, as in 
the case of the NN interaction in the '^Si~'^ Di channel, 
limiting how far the transformation can be run. When 
the SRG cutoff A — 1/si approaches the momentum 
scale at which a bound-state emerges, the strength of 
the SRG evolved interactions increases significantly. This 
happens because the transformation tends to move the 
bound-state eigenvalue to the low-momentum part of the 
Hamiltonian's diagonal, forcing the interaction to grow in 
order to maintain the bound-state at the right value. In 
the appendix we show that actually in the limit A — )■ 
the SRG equation becomes stationary when the SRG po- 
tential becomes the reaction matrix which in the single- 
channel case becomes 

limF.(p,p) = -^^, (36) 

A-i-O p 

and a similar equation is obtained in the coupled-channel 
case. Note that this diverges when the phase shift goes 
through 90°, a situation which only occurs for the 
channel for p ^ 100 MeV. On the other hand, there 



is no divergence problem when using Wegner's genera- 
tor T] — [Ds,Hs], because the bound-state eigenvalue is 
kept at the natural momentum scale as the SRG cutoff 
is lowered suggesting that actually the generator initially 
proposed by Wegner has better infrared properties. In 
this sense it would be interesting to check the present re- 
sults for the Wegner flow case. The advantages of other 
SRG generators have been considered recently '47 . 

Finally, we can fine tune the SRG cutoff so that we 
obtain the best possible fulfillment of the Wigner sym- 
metry, which is slightly above A. We call this Avvigncr- 
This, of course, would have far reaching consequences for 
the analysis of finite nuclei on the basis of symmetry. 

Our results for the 5-waves are presented in Fig. H] 
where the similarity cutoff runs from 5 fm~^ to 2 fm~ . 
We can see that both S'-waves evolve in the same direc- 
tion, becoming more attractive. At Awigncr = 3 fm~^, 
the SRG-evolved interaction in both S'-waves are practi- 
cally identical and for A = 2 fm~^ the evolved potential 
in the triplet channel is stronger than in the singlet case. 
As expected, the evolution in triplet state is faster than 
in the singlet state. The deuteron pole is at an imagi- 
nary momentum scale pd ~ i 0.23 fm~^ while the pole 
corresponding to the ^Sq virtual bound-state is at a much 
smaller momentum scale py ^ i 0.04 fm~l. Thus, when 
evolving to similarity cutoffs in the range we are consid- 
ering, the enhancement in the strength of the interaction 
that comes from using the simple generator 77 — [T, H{s)] 
is sensible only in the ^Si channel. As we see, at Awigncr 
the agreement between the and ^5'o SRG-evolved in- 
teractions is indeed remarkable both for the diagonal and 
the fully off-diagonal elements. It is important to note 
that the difference between the ^Sq and the '^Si phase- 
shifts evaluated from the initial unevolved potentials is 
reproduced at any rate through the SRG evolution, since 
the unitary SRG transformation changes the interactions 
(independently in each partial-wave channel) while pre- 
serving the corresponding phase-shifts. 

For the ZJ-waves the results are displayed in Fig. 
where we can see that generally for p',p < 0.7 fm^ 
Wigner symmetry is well fulfilled. Note that the better 
fulfillment of the Wigner symmetry for slightly large val- 
ues of p occurs for a similarity cutoff of 2 fm^^. A some- 
what similar situation occurs also for G-waves shown in 
Fig. [6j where now a better fulfillment of the Wigner sym- 
metry for a similarity cutoff of 4 fm~^. It is interesting 
to note that the value of the similarity cutoff where we 
have the better fulfillment of the Wigner symmetry is not 
unique, being different in each of the even angular mo- 
mentum waves we considered. Of course, while we do not 
expect a priori an exact fulfillment of 5t/(4)-spin-isospin 
symmetry it is remarkable how well it works taking into 
account that high-quality potentials have not been de- 
signed to implement the symmetry by any means. 

With the definitions of the potential we get the simple 
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expressions 

y-Lc = '^^J dzPL(z) [Vc - 3Vs + t{Wc ~ 3Ws)] , 

(37) 

V.Lc = 27ry dzPLiz) [Vc + Vs + t{Wc + Ws)] , 

(38) 



The non-triviality of the Serber symmetry is realized 
by noting that we just change the angle without changing 
the nucleons. To elaborate a bit further let us consider 
the standard Fierz identities. For the spin-isospin oper- 



where t = 4/ — 3. 
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ators one gets 



Pi 



1'2' 



-Pl'2'(0"l + 0-2) 
-Pl'2'(0"l + 0-2) 



[1- 



2 [1 - 3.7] , 

-((Tl + 0-2) , 
-((Tl + 0-2) , 



(39) 



and similar equations for isospin operators. These iden- 
ties imply that the following combinations of scalar func- 
tions are respectively even or odd in z 

[Vc - iVs + t{Wc - iWs)] 
= [Vc - iVs + t{Wc - iWs)] 
[Vc + Vs + t{Wc + Ws)] 

— Z 

= -[Vc + Vs + t{Wc + Ws)] (40) 

z 

Therefore, the orbital parity of the integrand and the 
Legendre polynomial in Eqs. (37) and (38) are the same 



and do not suggest that any of them vanish. Instead, the 
symmetries imply 

= Vc + Vs + Wc + Ws 
Vc-iVs + Wc-iWs = Vc + Vs-iWc-iWs 

(41) 

We remind that the large iV^ analysis of Ref. [53] yields 
Vc,Ws = 0{Nc) whereas Wc, Vs = 0{N-'^), suggesting 
from QCD that Wigner symmetry holds (second line of 
Eq. (41) in the even L partial waves, exactly as we do 
here for the SRG model. Note also that Serber symme- 
try under those conditions implies Vc + Ws — 0. These 
conclusions have been analyzed in great detail within sev- 
eral viewpoints [2SH2H] and are confirmed here within the 
SRG. In any case, our results suggest that it must be pos- 
sible to impose these approximate symmetries to a 
potential from the very beginning. 



V. CONCLUSIONS AND OUTLOOK 

In the present paper we have analyzed the concept 
of long distance symmetries as applied to the Similar- 
ity Renormalization Group. We have shown that very 
similarly to the case of the Viowk formulation, the sym- 
metry pattern of Wigner symmetry in partial waves with 
even angular momentum as well as the Serber symme- 
try in odd partial waves holds for the SRG cutoff around 
A ^ 3 fm^^. This is somewhat remarkable since SRG 
only provides an exponential decoupling between low- 
and high-energy modes bringing the effective interaction 



to the diagonal form, whereas T4owk corresponds to inte- 
grate out high-energy components. It is also notewor- 
thy that the symmetry arises at scales where in few- 
body calculations the induced three-body forces become 
small [Ulin]. 

In this work we have not evolved the three-body force 
through the SRG, so that no statements on the three- 
body SU(4) violating terms can be made, even know- 
ing that in the two-body case the SU(4) violating terms 
get small. From this point of view the extension of the 
present results to three- and four-body forces and the 
analysis of their symmetry structure would be of great 
interest. 

The underlying symmetry pattern unveiled by our 
SRG analysis appears intriguing and unexpected from 
the modern viewpoint of coarse graining high-quality in- 
teractions. From such a perspective this fits somewhat 
the vague concept of an accidental symmetry. Note that 
here we use the standard concept of accidental symmetry 
from quantum mechanics instead of the field theory one: 
a symmetry which is realized although not foreseen. 

There is a long tradition on the phenomenological con- 
sequences of Wigner symmetry in the properties of nuclei 
and nuclear matter (for a review see e.g. Ref. @H]). A 
recent work [49] analyses the SU(4) pattern of pairing 
forces within a Viowk framework, which quite naturally 
follow the symmetry pattern. Our results, in particu- 
lar the existence of an SRG scale at which the Wigner 
symmetry becomes quite accurate, not only provide a 
natural explanation for this fact but suggests to pursue 
the study further in future work within the current SRG 
framework. 

From a fundamental viewpoint, QCD large Nc based 
arguments foresee fulfilling the symmetry with a relative 
OiXjNl) accuracy, so one does not expect a perfect ful- 
fillment of the Wigner symmetry. On the other hand, 
nowhere in the design and optimization of the modern 
high-quality interactions which have been successfully 
applied to the structure of finite nuclei was the Wigner 
symmetry pattern explicitly implemented. In this regard, 
the accuracy with which by choosing a suitable SRG scale 
the symmetry seems to hold suggests that this is a prop- 
erty of the data themselves which emerges when the in- 
teraction is resolved with a specific length scale and not 
so much on the original (bare) potentials used to fit the 
A'^A'' scattering database. Finally, one should recognize 
that while the use of symmetries for coarse grained ef- 
fective interactions is not mandatory we expect by ex- 
plicit symmetry considerations additional simplifications 
of the complicated nuclear many-body problem. Work 
along these lines is in progress. 
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so that the potential may be written as 

V = Vc + (2C|^(4) - l?>)Ws . (AlO) 

The particular choice Wg = corresponds to a spin- 

isospin independent potential, but in this case no dis- 
tinction between the 6^ and lOg supermultiplets arises. 



Appendix A: Wigner symmetry for NN 

For completeness we remind here some features of the 
Wigner SU{A) spin-isospin symmetry. It consists of the 
following 15-generators 



rpa \ ^ 

A 

^ A 

= X 2^(^at^a^ 



(Al) 
(A2) 
(A3) 



where and (t\ are isospin and spin Pauli matrices for 
nuclcon A respectively, and is the total isospin, 5* the 
total spin and the Gamow- Teller transition operator. 
The quadratic Casimir operator reads 



r~i rparp I C I /^ia/^ 



(A4) 



and a complete set of commuting operators can be taken 
to be Cs[/(4), T3 and Sz,Gz3- The fundamental repre- 
sentation has Csu{4) = 4 and corresponds to a single 
nucleon state with a quartet of states p ^, p I, n ^, n ^, 
with total spin 5* = 1/2 and isospin T = 1/2 represented 
4 = (5,r) = (1/2,1/2). For two-nucleon states with 
good spin S and good isospin T Pauli principle requires 
(^—"jS+T+L _ -with L the angular momentum, thus 



Csu(4) = ^(a + r + ar) + 



15 

y 



(A5) 



where t = ti ■ T2 = 2T{T + 1) — 3 and a = ai ■ a2 = 
2S{S + l)-3. 

One has two SU(4) supermultiplets, which Casimir val- 



ues are 



^SU(4) ~ ^SU{4,) — ^ ' 

/^01 _ /^lO _ r 

*-^Si7(4) — *-^Si7(4) — <^ ' 



(A6) 
(A7) 



corresponding to an antisymmetric sextet 6a = (0, 1) ® 
(1,0) when L = even and a symmetric decuplet IO5 = 
(0, 0) ® (1, 1) when L =odd. According to LSJ quantum 
numbers we have the following supermultiplets 

(^50,^ 5i) , (lPi,3 Po,l,2) , ('^^2," i5l,2,3) . . . (A8) 

When applied to the; NN potential, the requirement of 
Wigner symmetry for all states, implies 

Vt = Wt = Vls = Wls = 0, 
Ws = Vs = Wc, (A9) 



Appendix B: The infrared cut-off limit of the SRG 

In this appendix we show that when the kinetic en- 
ergy is taken as the generator of the SRG transforma- 
tions the SRG evolved potential becomes the standard 
reaction matrix, i.e. 



(Bl) 



when there arc no bound-states. To avoid unnecessary 
mathematical complications, we will analyze the problem 
in the discretized form like the gauss integration points 
used in the numerical calculation. The SRG equation in 
the basis where T is diagonal becomes 

^ = -{Ei - EkfVik + + Ek- 2Ei)ViiViiB2) 



The discrete representation has many advantages. One 
can see that along evolution in the discrete representa- 
tion one has an infinite number of constants of motion 
j^Ti{V") = due to the commutator structure of the 
SRG equation ^. 

The fixed-point solution implies ^^j^ which requires 
that [[T, V],V]=Oso that V and [T, V] become diagonal 
in the same basis, not necessarily the one where T is 
diagonal, i.e. 



Va^ = Vada0 (B3) 

= ^(T,^F^^-F„^T^;3)=T„^(F^-14)(B4) 

7 

The second line requires that for a ^ /3 then Tq,^ = 
which means that T is diagonal also. Therefore Vij = 
ViSij , i.e. V is diagonal in the basis where T is also 
diagonal. If we write now the Lippmann-Schwinger (LS) 
equation for the reaction matrix 



R{p',p;k) = Vip',p) 



+ 



2 J° 

71" 7o 



fdq 



V{p',q)R{q,p;k) 



(B5) 



^ Mathematically, such a property is ill defined in the continuum 
limit since even for n = 1 one has Tr{Va) = /q°° p'^Va{p,p) = 
r'^drV{r,r) which for a local potential V{r,r') = V{r)S{r — 
r ) diverges as the momentum cutoff. Also, the trace of a com- 
mutator Tr[yl, , B] only vanishes when both Tr{AB) and Tr{BA) 
are finite, as the choice A = p and B = x clearly illustrates, since 
\p,x] = —ih and hence Tr[j),a;] = — i/lTr(l = 00. 
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FIG. 6: Comparison between diagonal, V{p,p), and fully off-diagonal, V{p,0), matrix-elements of the SRG-evolved potentials 
for the G- waves (in fm) as a function of the CM momentum p (in fm~^). 
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With this normahzation the phase-shift in the one- 
channel case reads 



tanS{p) 



P 



(B6) 



Note that the SRG actuaUy impUes that the phase-shift 
is constant along the evolution, so that one may take any 
Vs in Eq. (B5). The discrete version of this equation for 



the half off-shell i?-matrix, Rij — R{pi,pj;pi reads 

-^QQk^ 2 



(B7) 



where the principal value corresponds to exclude the 
point in the sum. Thus, for a diagonal potential we get 



(B8) 



Turning to the phase-shift and going to the continuum 
limit we thus obtain the assertion, so that Eq. (36) is 



obtained. The proof also holds for the continuum limit, 
provided a momentum cutoff in the integrals is supple- 
mented. 

The infrared solution suggests looking for perturba- 
tions around it. Actually we will see now under what 
conditions are these perturbations stable. If we write 

V,k{s) ^ VAj + AV,j (B9) 

so that we get to first order in the perturbation 

AV:,{s) = {E, ~ E,){E, H- - E, - \/,)Ay,,(s)(B10) 

which yields the solution 



(Bll) 



The diagonal part is constant as required by the prop- 
erty Try(s) = const.. To identify this contribution we 
use again the LS equation and find that /SVa — and 
AT/ij(oo) = Rij for i ^ j. Thus, we get 



Al/,,(s) = (l-%)i?,,e 



-s{E,-Ej)(E,+V,-Ej--Vj) 



(B12) 



which clearly shows that the fixed-point is stable pro- 
vided {Ei — Ej)Vi — Vj) > 0. Departures from it measure 
some off-shellness of the potential. Going to the contin- 
uum limit we get for p' ^ p 



where for large p we obtain R{p,p) — > —S{p)/p = V{p,p) 
which means that at high-momentum the kinetic energy 
dominates and the fixed-point is stable. If there are 
no bound-states, we have that 6{p) never becomes 7r/2 
(Levinson's theorem states that 5{0) — S{oo) = nsTr, with 
n B the number of bound-states) . However, for a pole at 
p = Po we get 



R{p,p) 



1 



reg. 



(B14) 



Po{p - Po)5'{pg) 
where reg. means regular contributions, so that for p > 

Po > p' and 5'{po) > or p' > po > p and S'{pq) < 
becomes stable and unstable otherwise. In the '^Si chan- 
nel one has 5'{po) < which means that the corrections 
increase dramatically for p > po > p' . 
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